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Abstract 

Bound state problem in the relativistic QED is investigated by the 
functional integral methods. The ortho- para- positron mass differ- 
ence is calculated. Contribution of the " nonphysical" time variable 
turned out to be important and leads to the nonanalytic dependence 
of the bound state mass of the order as. It is shown that the relativis- 
tic and non-relativistic QED gives different results for this mass shift. 
In addition so-called abnormal states as "time excitations" arise. Se- 
quential application of relativistic QED to bound state problem is in 
contradiction with real ortho- and para- positronium bound states. 

The conclusion: the relativistic QED is not suited to describe real 
bound states correctly. 
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1 Introduction. 



We believe that the relativistic quantum electrodynamics (QED) is a uniquely 
correct universal theory giving an exhaustive description of all interactions 
between electrons and photons including possible bound states like positro- 
nium. Only our inability to calculate something out of perturbation method 
does not permit us to obtain all the desired details. Earlier, some scientists 
considered that QED should have its own applicability region. A short review 
of the history and the development of quantum field theory is done in pQ. 
Supporting these doubts we will show in this paper that the sequential use 
of the standard QED does not give a correct description of the positronium 
spectra, namely, the ortho- para- positronium mass difference. 

First of all let us realize what is the status of bound states in non- 
relativistic quantum mechanics (QM) and relativistic quantum field theory 
(QFT). In what follows, we restrict ourselves to the discussion of positronium 
in QED. The difference between QM and QFT is shown in Table I. Let us 
give some comments. 

The total Hamiltonian H = Ho+gHi can be constructed in QM and QFT. 
However, in QM H is a well defined operator, so that the non-relativistic 
Schrodinger equation is mathematically correct and time development of a 
quantum system can be described. Solutions of the Schrodinger equation 
contain both free and bound states. One can remark that in QM a bound 
state (positronium) is created by real particles (electron-positron), i.e. con- 
stituent particles are on mass shell but are no virtual particles. 

In QFT the Fock space T is defined by the noninteracting free Hamilto- 
nian H Q and contains the free particles only. However, Hj is not defined on 
J 7 . As a result, the bound state as an eigenvalue problem of the relativistic 
Schrodinger equation on the Fock space cannot be formulated mathemati- 
cally in a correct way (see [2]). Besides, the time development of quantum 
field system cannot be obtained. The only way to overcome these problems 
is to construct the S-matrix which contains all elastic and inelastic scattering 
amplitudes of free particles from the time t — > — oo to the time t — > oo. It 
is important that the S'-matrix is a unitary operator on the Fock space. It 
means that the bound states like positronium, which is a unstable particle, 
cannot belong to any Fock space in principle. In addition, our computing 
abilities are restricted to the perturbation theory. 

Nevertheless, we believe that the S'-matrix amplitudes should contain 
some correct information on possible bound states. The simplest way to re- 
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alize this idea is to postulate that a bound state is a simple pole of an elastic 
scattering amplitude of constituent particles with appropriate quantum num- 
bers. It means, that the amplitudes out of mass shell and out of perturbation 
approach should be calculated. Standard methods to go out of perturbation 
calculations are reduced to sum appropriate classes of Feynman diagrams 
and this summation can be formulated in a form of integral equations. The 
best known approaches are the Bethe-Salpeter and Schwinger- Dyson equa- 
tions. There is numerous literature devoted these equations (see, for example, 
[3J HI El El [7J HI E]). The important difference comparable with the nonrela- 
tivistic case is that bound states in these equations are created by particles 
which are out of mass shell so that the role of time becomes important. 

One remark on these equations. We know that the perturbation series 
are asymptotic series so that the problem is how to sum them? The exact 
amplitudes should have some singularity at the point a = in QED (see 
[10J). What is a precise character of this singularity is not known up to now. 
Standard perturbation expansions are connected with Feynman diagrams. 
Usual methods are reduced to summation of an appropriate class of Feynman 
diagrams. Result of a summation of any definite class of Feynman diagrams is 
a kind of geometrical progression, i.e. it is an analytic function at the point 
a = 0. However, it should be stressed that the generally accepted point 
of view - non-perturbed behavior is a sum of a definite class of Feynman 
diagrams - is not true. 

One of probably successful proposals to calculate the relativistic correc- 
tions to bound state problem is the so-called non-relativistic QED (NRQED) 
(see [H]). The basic idea is that the QM is correct, only non-relativistic 
momenta are responsible for bound state properties. In other words, the 
Hamiltonian should not depend on time and the problem is to find somehow 
relativistically small corrections to the non-relativistic Coulomb potential. 
The basic idea is that for small coupling constants the Born approximation 
is a good approximation which is directly defined by the Fourier transform of 
the potential. The aim is to extract from the relativistic S-matrix some rel- 
ativistic corrections to non-relativistic Hamiltonian. The hypothesis is that 
the scattering amplitudes in the non-relativistic Schrodinger theory and the 
relativistic S'-matrix theory should coincide in the low energy limit. The 
procedure is to write down the non-relativistic Lagrangian with a set of all 
possible terms, and coefficients in front of them are calculated by identifica- 
tion with appropriate amplitudes of relativistic S'-matrix. This prescription 
allows one to remove effectively time out of the relativistic equations, in 
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other words, to place all intermediate particles on their mass shell. It seems 
NRQED is supported by experimental data. 

Another quantum field idea is that a bound state is defined by an asymp- 
totic behavior of the vacuum mean value of the corresponding relativistic 
currents (see, for example, [12]) with desired quantum numbers: 

(0|J(x)J(0)|0) = £ (°I J 0*> HJ(0)|0) = J2e- En]xl \ (0|J(0)|n) | 2 

n n 

~ e~ Mmin ^\ (0|J(0)|mm) | 2 for \x\ -»> oo. (1) 

This formula gives a possibility to calculate the mass of the lowest bound 
state | mm) if M min < 2m. Essentially, the space of states {\n)} is supposed 
to contain possible bound states although we saw that the Fock space can- 
not contain unstable bound states. These vacuum mean values (EJ) can be 
represented in closed forms by functional methods. The functional meth- 
ods permit one to get formally the exact representations for Green functions 
which are not connected directly with Feynman diagrams, so that it is possi- 
ble to go out of standard perturbation expansions using asymptotic methods. 
Development of functional methods permits one to get the exact character of 
non-analyticity at the point a = and to clarify the role of "time" in bound 
state formation. Exactly this approach will be used in this paper. 

The practically unique experimental object to investigate the bound state 
problems is the positronium which is the result of pure QED interaction. On 
the one hand, the positronium is not a stable state. It cannot belong to 
the asymptotic Fock space. Nevertheless, it exists. The binding energy of 
positronium itself is not measured with great accuracy but the mass difference 
of two possible states, ortho-positronium (l 3 5i) and para-positronium (l 1 ^), 
is known with very large accuracy 

Ae = - e para = 203.38910 GHz = 8.4115 • 10" 4 eV (2) 
7 

= 0.580487 a 4 m e = — a 4 m e • 0.99512... 
12 

The main contribution can be explained by the non-relativistic Breit poten- 
tial approach (see, for example, [321 EU EI3 EE] ) taking into account scattering 
and annihilation channels 
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If we apply the relativistic current formula ([T]) to the positronium prob- 
lem, we can write 

(0|J(x)J(0)|0) = ]T e- lE "^\ (0|J(0)|n)| 2 + ]T e ~^\\ (0|J(0)|n) | 2 

particles photons 

where 

e~ iEnlxl \ (0|J(0)|n) | 2 ~ e " M " N 

particles 

and annihilation channel looks like 

1 



Y e- iB " N |(0|J(0)ln>| 2 

\x\ 



photons 

It means that the annihilation channel does not take part in the bound state 
formation in contradiction with the non-relativistic potential approach. 

Another point: we want to understand what is the role of TIME in for- 
mation of bound states. 

In this paper we apply functional methods to calculate the asymptotic 
behavior of vacuum mean value ([1]) of relativistic currents for positronium 
and clarify the role of time in the formation of bound states. 



2 Lagrangian and bound states 



All our calculations will be performed in the Euclidean space. The La- 
grangian of the electron field if) and the electromagnetic photon field 
looks like 

1 



(3) 



dpAvix) - d v A ll (x). 
The electron propagator has the standard form 



S(x — x) = (T 
The propagator of the photon vector field is 

D^{x -y) = (A^x)A u (y)) = S^ u D(x - y) + 



dp e 



ip(x— x') 



(4) 



(27r) 4 m — ip 

D d (x-y), (5) 



2 



D(x) 



dk e lkx 1 



(2vr) 4 k 2 {2n) 2 x 2 ' 



DJx) 



dk 



d{k 2 



(27T) 



k 2 k 2 
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2.1 Two-point Green function 

The object of our interest is the gauge invariant two-point Green function 



G 



-m]ip) 



■ @(x)ri>(x))@(y)rr(>(y)) (6) 

Here T is a Dirac matrix which defines the local vertex with quantum numbers 
of the state Jr = (ipTip). We have for para-positronium T = ij 5 and for 
ortho-positronium T = 7^. 

After integration over the electron fields ip and ip we get 



where 



Br(x 



Gr(x - y) = B r (x - y) + U r (x - y), 



y ) = J ^ e -j(^^)+T[A] .Tr[rS(x,y\A)rS(y,x\A)\, (8) 



(7) 



and 



H, 



ix-y) = J 



C 



T![TS{x,x\A)\.Tr[rS{y,y\A)\. 



Here S(x, y\A) is the electron propagator in the external field A^. 

1 



S(x,y\A) 



i(p + eA(x)) — m 



5(x - y) 



(9) 



The functional 



T[A] = Tr In - — ^ = Tr In 



ip — m 

A 1 A 1 

A— A— — 



1 + ieA- 



1 



ip — m 



ip — m ip — m 



+ 0(eM 4 ) 



e 



S (x - y) 



|| A„(x)Il^(x - y)A v {y) + 0(e 4 A 4 ), 

dp e - ip{x - y) 



1 

ip — m 



-S(x - y) 



(27r) 4 ip — m 



IV(x - y) = Tr [7^0(2; - y)luS Q (y - a;)] , 
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13 ( x - y) ~ e -M|x-y| 




Figure 1: Terms B and H 

describes radiation corrections to the photon propagator and to the photon- 
photon interaction. In this paper, we neglect this term because it does not 
contain spin-spin interaction and, therefore, does not contribute to ortho- 
para- positronium mass difference in the lowest corrections. 

The loop Bp contains all possible (?/T?/))-bound states. If the mass of the 
lowest state Mr < 2m, then the asymptotic behavior of this loop for large 
\x — y\ looks like 

B r (x - y) ~ e ~ M ^-y\ (10) 

where Mr is the mass of the lowest state in the current (ipTip), i.e. the mass 
of a possible bound state. This mass can be calculated by the formula 

M r = - lim -!- In B r (x) = 2m - e r . (11) 

\X\ 

Here er defines the binding energy of the lowest bound state. Our aim is to 
calculate the functional integral (IE]) in the limit \x — y\ — > oo and to find Mr, 
according to ( TTTT) . 
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The loop H describes so the called annihilation channel and contains 
long-range contributions of photons: 



H r (a; - y) ~ - ^ 

\x — y\ 2 

This term does not contain any bound state. 

Graphic representations of the loops Bp and Hp are shown on Fig.l. 

3 The electron propagator 

The propagator of the electron fermion field satisfies the equation 

[i(p + eA(x)) - m]S(x, y\A) = 5(x - y), (12) 
For the gauge transformation 

Ap{x) — ► A,(x) + dJ(x) 

it is transformed as 

S{x, y\A + df) = e iefix) S{x, y\A)e- ief(v \ 

so that the loop (jSJ) is gauge invariant. 

The solution of the equation (fl2|) can be represented by the functional 
integral (see, for example, [17]): 

1 

S(x, y\A) = - 5(x — y) 

i(p + eA(x)) — m 

= \i(Px + eA(x)) + m] ■ —, r-r- — — -5(x — y), 

lU y " J (p + eA(x)y + ^F^x) + m 2 1 Uh 

00 r %2 

= [2(p x + eA(a;)) + m\ J 



s 2 s 

D?7 - /(ft^+ie /ii!i M (t)jl„(«(t)) I f f ' dta^(t)F^(z{t)) 

'e T ( e k (13) 



C 



The boundary conditions are r/(0) = rj(a) = and the normalization is 

Drj \ f , fj 2 (t) 



The symbol Tt means the time-ordering of the matrix cr M!/ (r) to the time 
variable £. 

The representation (fl3|) is quite close to functional representation of the 
propagator for a scalar charged particle (see [9]). The main functional struc- 
ture is the same. 

The representation (Tl3l) is obviously gauge covariant because 



d 

5 J dt ^(t)A^z(t)) = J dt z^t) — f(z(t)) 

^ 

dt^f(z(t)) = f(z(s)) - f(z(0)) = f(x) - f(y). 
o 

As it was said above, our aim is to calculate the functional integral (jSJ) 
in the limit \x\ — > oo (we put y — 0). We want to calculate the main 
contributions to the binding energy assuming the coupling constant a to be 
small. In this case, for large x — > oo and small a the saddle-point in the 
integral over s is realized for s — — . Putting; 

x = (x, x 4 ) =>■ (0, x 4 ), =>- X4, = X > 0, t = — r, 



??? 



one can get for X — > oo 



S(x, 0\A) ^(1 + 7o)e" mX • 5(x) 

5(0, =► — ^(1 - 7o )e- mX ■ S(x) (14) 

A 2 

R[z) =T T le ° 



with 

z{T) = nT + r]{T) = S (is) 
[ r + rji{T). 

We shall use this representation in what follows. 
3.1 Mass of the bound state 

The next step is to substitute electron propagators S(x,0\A) and S(0,x\A) 
in the form (THI) into the representation ([8]) for the Green function Br(x) and 
then to integrate over the photon field A. We have for large X — > oo 

X 

B T (X)~e- 2 ™ x ll^^e o JtfX.mM (16) 

with 

(17) 

DA / rfr4 1) (r)A M ( 2 ( 1 )(r))+ie / dri^ (r)A M ( 2 ( 2 > (r)) 
p 

c 

r (1 + 7o ) J2[z«] r (1 - 7 o) R[z {2) ] 



The mass M-p of the bound state with quantum number T is defined by 
the formula (flip . 

The integral (|17p over the photon field A can be calculated explicitly. 
The result of the calculations is shown in Fig. 2. We will not write down 
this simple long formula. Our aim is to find the ortho-para mass difference 
in the lowest approximation of the functional method. Therefore, we omit 
all terms connected with contributions to the electron propagator in the 
loop and take into account the dominant terms responsible for positronium 
formation and desired ortho-para mass difference (see Fig. 3). We get in the 
lowest approximation over spin-spin interaction 



7r[X; 771,772] 



DA 



iejdr 4 1) (r)A fl ( Z ( 1 )(r))+ie]'dr (r)A^(z^ (r)) 
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Figure 2: All diagrams contributing to (IT 



■i 



iTr|r(l + 7 o)r(l-7o) 



.2 V 



+ IoW JJ dTldT2 ' r(7 ° + ^^rC-To + • i^(z (1 Hn))iV(* (2) fa)) 



Here 

4 0) = ^Trr(i + 7o )r(i- 7o ). 

The main functional responsible for the bound state formation looks like 

x 

W[X; VuV2 ] = e 2 JJ d n dr 2 ^\r^{r^)D^\n) - *«(r 2 )) 
o 

X 

e 2 [[ dndr 2 ^\r x )zf{^)D{z^{r x ) - z^{r 2 )) (18) 



' o 
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Figure 3: Diagrams which are responsible for the bound state 

dk e ifc(*W(n)-*«(7a)) 



,2 



dndr 2 e\ n )e\r 2 ) 



(2tt) 4 k 2 
The functional responsible for the ortho- para- mass difference is 



2 X I 



K r [X; r/!, r/ 2 ] = ^— ^ J J dr x dT 2 ■ -Tr T(l + 7„)<7 w r(l - -y )a up 

o 

°~ D(zW(n)-^(r 2 )) (19) 



o v ' 

where 

£ r (fc) = ^ . ^Tr r(l + 7 o)<V,r(l - l0 )a pv = 4°> ■ A T (k). (20) 

The term W is responsible for positronium bound states. The term 
describes the spin-spin interaction and it is responsible for the ortho- and 

2 

para- mass difference. All neglected terms give the next to a — |- perturba- 
tion contributions. 

Let us introduce the notation 

X 

w DrfrDrh ~f / dt[m 2 (t)+f,l(t)]+w[x m , V2 ] 
~~C 



d< V2 = ^p^e o (21 ) 
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with 

J (X) = J dal m = Je^ x (22) 

then one can write 

B r (X) = e- 2mX J I dal m ■ {4 0) + K r [X; % , m ] + 0(e 4 )} 

= J £(?) e - (2m - e °- £r)x (23) 

The binding energy eo does not depend on the spin of electron-positron sys- 
tem and it is defined mainly by the Coulomb interaction. The binding energy 
er depends on the spin of constituents and defines the ortho-para mass dif- 
ference. It looks like 

e 2 1 rr , . r dk / oik{zW(Ti) _ z (2) (T2)) 



e r = lim -— // dndr 2 / — — • (e^* v [Tl) ~ z (T2 >>) A r (k), 

x^oo Am 2 XJJ J 2vr 4 \ 'mm y h 

o v ' 

(24) 

where the average is 

\ 6 I mm ~ MX) 11 mm [ ] 

The functional W in (|21[) contains terms defining the Coulomb and "time" 
interactions. The measure da^ V2 contains the space r\ and "time" 774 func- 
tional variables. Our aim is to evaluate the contribution of "time" interaction 
to the bound state formation. Thus our direct problem is to calculate the 
integral (T25D . 



4 Para- and ortho-positronium mass differ- 
ence 

In the quantum field theory there is the problems how to define the bound 
state. The point is that the relativistic currents of the type (ipTi/j) are classi- 
fied by the relativistic group in the space R 4 while the physical bound states 
are classified by the non- relativistic group in the space R 3 . As a result the 
physical states are described by an appropriate mixture of different compo- 
nents of different relativistic currents. In addition, the angles of mixture 
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are not known a priori. In Table 2 the non-relativistic quantum numbers of 
different components of relativistic currents are listed. 

One can see that each current (ipTip) with quantum number J p is deter- 
mined by a mixture of two relativistic currents 

S(0 + ): r 5 = / cos S + 7o sin S , 

A(l + ) : T A = 757COS0A + i[7 x 7]sin0 A , 

V(l~) : IV = 7cos0y + Z7o7sin0y, 

P(0~) : T P = ^75 cos0p + 7 5 7o sin0p. (26) 

Generally speaking,in order to define the angles 0s, a, Oy, Op, some addi- 
tional argumentation should be used. However, we shall see that the desired 
masses do not depend on these angles at least in the lowest approximation. 

Let us come to formula (1201) . The results for and Sr(^) with currents 
(|2"6"|) are listed in Table 3. One can see that the states S^O" 1 ") and A(l + ) do 
not exist at all. The masses of bound states P{0~) and V(l~) do not depend 
on the mixing angles 6p and Oy. 

According to fT2"4"j) . the desired mass difference is defined by the formula 

5M = e v - € P 

v ' 

= - ■ - lim - // dr ld r 2 [[ ■ e ^, 2) V (2?) 

3 4 x^oo X J J 1 2 JJ (2tt) 4 k 2 + k\ K ' 

5 The lowest contribution 

The lowest main contribution to positronium bound state is defined by the 
integrals f l2"2"j) and f l23|) . It is convenient to extract the a 2 dependence in order 
to extract the non-relativistic Coulomb potential (see [9]). For this aim in 
the representation fZTj) let us introduce the new variables: 

Y = a 2 mX, r = — - — , = a 2 mq, k = amq. 
a 2 m 

v(t) = —${v), ri(t) = — 

am am 
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The parameters X and Y are infinitely large quantities. 
We get 



Jo(Y) 



— e 

(28 



C 



with 



Y 



= J J dv x dv 2 (l + aCi(^i)) (l + a&faa)) + a 2 £iK)£ 2 ( v 2) 
o 

// 47r 3 q 2 + a 2 g 2 



(29) 



It is important to stress that the functional £ 2 , £i, £2; a ] is n °t ana- 

lytic at the point a = 0, so that the relativistic corrections cannot be obtained 
by a regular method. Calculation of the integral ( 128|) is not a simple prob- 
lem. Therefore, we restrict ourselves to calculation of the next relativistic 
correction to a, so that we neglect terms with a in square brackets in (I2"9~j) 
and introduce the variables 

&(r) = R(r) + ~p(r), £ 2 (r) = R(r) - \p{r), 

ei(r) = i2(r) + ip(r), £ 2 (r) = i2(r) - ^(r), 

The variables R, p and i?, p describe the center of mass and relative coor- 
dinates in configuration and "time" spaces, respectively. We shall see that 
the "time" variables give the important contribution to the desired mass 
correction. We get 



Jo(Y) 



J J DpDRDpDR -J dv [kP (v)+R, 2 (v)+^p 2 {v)+R 2 (v)]+W[R„p,R,p;c 



C 



where 



1 

W[R, p, R, p; a] = J J dvidv 2 



(30) 



(31) 
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// 



dqdq e i 9(^-^+a(i?(i'i)-i?(^2) + i(p(i'i)+p(«2)))+iq(R(t'i)-R(^2) + |(PK)+P(t'2)) 

47T 3 q 2 + a 2 q 2 



The variables R and R are connected with the continuous spectrum so 
that in the lowest perturbation order over the interaction functional W[R, p, R, p; a] 
we can integrate over variables R and R. We get for large X: 



MY) 



DpDp -]dv[\?{v)+\p 2 (v)}+W[p,p-M 
p 

C 



(32) 



Here 



W[p,p; a] = jj 
dv\dv2 II 



DRDR -!*>[&<?)+#<?)] 
e 

c 



W[R, p, R,p; a] 



JnJn e,i<l(vi-V2)-2[a 2 q 2 +<i 2 ]\vi-V2\ 

d ^ d< i e 4 . e -^(p(vi)+ P (v2))-^ (P(«l)+P(«2))_ 

47T 3 q 2 + a 2 g 2 



The functional W [p, p; a] is not analytic at the point a = 0. Nevertheless, 
one can extract the lowest terms to a. 



W [p, p; a] W [p, p; a] 

Y 

= II dv\dv2 



dcjdq ^ qfa _„ 2 )-i q 2 kl _„ 2 | _ e _«a Wt)1 ) +p ( W2 ))_^(p( t;i ) + p( W2 )) 
47r 3 q 2 



/■ dq 




/ 27r 2 q 2 7 l 












/ 2tt 2 


7 q^ 
.0 



(if 



1 - f p'(^) 



-^\p(v)\-idP(v) 



Y 



Taking into account the correlations (see Appendix III) one can get 



5(p(v)) = (5(p(v))) p + [5(p(v)) - (8(p(v))) p 



8tt 



+ 



5( P (v)) - 



1 ■ 

8tt. 
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because 



Finally, we have 



(S(P(V))) P = ^(0) 



Y A Y 

av a 



8tt 



\p(v 

u 

The measure in fl32l) is reduced to 



Here 



da 



DpDp -]dv[\f(v)+\p\v)]+W[p,p-M 



C 



-e o 



— t- da = dapda p , 



da 



Dp -S dv [kP 
P 



P = ^e o 



(34) 



Dp -J*>[*?W+fW)\] 
da p = —e ° 

This measure consists of two components - the standard non-relativistic 
Coulomb potential term and the "time" term which is not taken into account 
in any usual calculations. This term corresponds to one-dimension linear 
potential with coupling constant a and leads to non- analytical behavior of 
energy on a. 

One can calculate (see Appendix III) 



J, 



Coulomb 



da, 



Dp -J dv 



£(2£ + l)e- 



-X- 



C, 



l^(0)| 



Jtime = da p = -fe - = E e Se "m|$„(0)|' 



(35) 



17 



Thus the general spectrum of the Coulomb and "time" potentials is 

„2 



EL 



a 



rn 



(36) 



As a result we have the "time excitations", or abnormal states, connected 
with the fourth component of 4-dimensional space. These states appear in 
solutions of the Bethe-Salpeter equation. Up to now it is not known exactly 
these states does or does not exist in reality. It is the second reason why 
relativistic QED does not describe correctly the real bound states. 

6 Mass difference 

The desired mass difference is defined by the formula (1271) which can be 
represented as 

8M= - ■ a 4 m ■ A(a), 
3 



Y 



A(a) = hm — / / dvidv- 
v ' y->oo Y J J 
o 



1 /■'/■ f [ d(\dq q 2 e i l( v l- v 2)~j[a 2 q 2 +q 2 ]\v 1 -v 2 \ 



+ a A q 



2 n 2 



-£(p(*i)+p(*a)) 



e -^r(p(«i)+p(«2)) 



p ^ i p 

The averaging over the fields p and p gives (see Appendix III) 



e 2 



lP(n)+p(r 2 ) 



e -*f (p(n)+p(r 2 )) 



K 

Finally we have 

A(a) = A I dr II dkdq ■ — — • e ^-^V+k 2 ]M (3?) 

7r d J J J k 2 + crg^ 



■l)^+l)C^(^]^(-l)"e * 



«,= 4— (e«-eo) 



00 n— 1 ? 

EEE(-i) £+K ^ K (a) 

n=l £=0 k=0 
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After integration over r and angles one can get 



A„i K (a) = —r 



7T" 



dkdq k* k 2 + a 2 q 2 + 1 - 4y + ah (e K - e ) 

J k 2 + a 2 q 2 + a 2 q 2 + 1 _J_ + a l {eK _ eo) y + 16q 

k 



(2£+l)C 



ni 



A K (azq) 



(38) 



The numerical results are shown in Table 4. For the function A (a) we get 

Ao(0) = EE(-l)'A^o(0) = 1.0, 

n=l e=o 

3 n-1 

Ao(a) = E E(- 1 )' A ^o(«) = 0.9641, 

n=l e=o 



3 n-1 2 

A(a) = E E E(- 1 )^ A ^(«) = 0.952754. 

71=1 1=0 K = 

Obviously, this result is in contradiction with the existing experimental 
number (A = 0.99512...). 

7 Breit potential approach 

One of the attempts to describe the bound state problem is the Breit potential 
approach (see, for example, p~3] ). Let us consider the elastic electron-positron 
scattering 

e p + e k =^> e p > + e k > 

The scattering amplitude in the lowest order of relativistic S'-matrix theory 
is described by the Feynman diagrams shown on Fig. 2 and looks like 

M = -e^uip'h^p^D^p-p'M-khM-k')] (39) 

+ e 2 [u(-k)^^u(p)]D^(p + k)[u(p)^ u u(-k')} 

where the first term (a) is connected with scattering and the second one (b) 
with annihilation channels. The spinors u(p') and u(k') are the solutions 
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(a) 



(b) 



Figure 4: Diagrams defining the Breit potential: (a) - scattering and (b) - 
annihilation channels. 

of the Dirac equation. This amplitude in the non-relativistic limit should 
coincide with the non-relativistic Born approximation which defines the ef- 
fective electron-positron potential. This potential should be introduced into 
the non-relativistic Schrodinger equation. 

The Dirac spinors u(p') and u(k') define the relativistic corrections to the 
non-relativistic potential. It should emphasize that electrons and positrons 
are on the mass shell, i.e. they are real physical particles. Thus the time is 
removed from the interaction Hamiltonian. 

The part of the Breit potential which is responsible for the ortho- and 
para- mass difference looks like 

U(r) = U sc (r) + U an (r) = — + -.—(*_CT + )S(r). (40) 

r 12 m z 

where 

UJr) = — + -•— (<r_<r + )tf(r), 

U an (r) = -•— (<T_<r+)<y(r). 

4 m z 

We want to stress that the coefficient ^ is the sum of the contributions from 
scattering and annihilation channels 

7 /1 
12 " 
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Taking into account 



T 9 . , am 

* "<°> = -ST- 

((<T_CT + )) ort foo — (( cr - cr +))para = 4, 

one can get for the ortho- and para- mass difference 

7 2air a 3 m 3 7 4 

^ £ ortho tpara TTT ' ^~ ' Z ' ^ TrT^ ^ 
12 TTT 87T 12 

This result is in good agreement with the experimental data and, therefore, 
supports the point of view that the annihilation channel plays the essential 
role in the formation of the positronium. 



8 Conclusion 

In conclusion, one can say that the functional approach is the best math- 
ematical representation to preserve the gauge invariance. The developed 
technique of calculations permits one to get accurate results in QED where 
the coupling constant a is small. The lowest approximation of this functional 
representation is the pure non-relativistic Feynman path integral representa- 
tion of the non-relativistic Schrodinger equation with the Coulomb potential. 
One can see that any regular series for next corrections to a do not exist and 
these corrections can not be reduced to some terms to the non-relativistic 
potential in the Schrodinger picture. In other words,the " nonphysical" time 
coordinate is important and leads to corrections which is not analytic of the 
order as. 

There exists a contradiction in the current algebra formula ([T]). On one 
hand, it is supposed that the space of states {|n)} can contain possible bound 
states. However on the other hand, in reality it is the Fock space of free 
electrons and photons which does not contain any unstable bound states. 
Nevertheless calculations of the functional representation for an appropriate 
Green function in the limit t — > oo indicate that a bound with Mb OU nd < 2m 
does exist really. Besides, the current algebra in QFT excludes influence of 
the annihilation channel for the bound state formation. 

Our calculations show that the role of time is very important and give 
essential contribution into bound state mass. The next radiation correc- 
tions, connected with time excitations, to electromagnetic mass difference to 

2 

positronium are of the order «3, i.e. they are to large. 
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In addition, the "time excitations", or abnormal states arise in QFT 
calculations but they are not exist in reality. 

The experimental value of ortho- para-positronium mass difference is de- 
scribed in the framework of the Breit potential picture with attraction of the 
annihilation channel. Thus, explanation of experimental value para- ortho- 
positronium mass difference requires to take into account annihilation chan- 
nel for effective potential. 

One can conclude that in the relativistic QED time corrections are im- 
portant, but the bound state problem requires the non-relativistic potential 
description where the time variable does not play any essential role. 

The conclusion: the relativistic QED is not suited to describe real bound 
states correctly. 
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9 Appendix I 

We use the following representation for 7 matrices: 

70= (J -l)' 7= (-l 0)' 75= (? I) 

and 

We should calculate the traces 

4 0) = ^Trr(i + 7o )r(i- 7o ) 

and 

Er ( fc ) = ^ • Itt r(l + 7 oK p r(l - lo )a up = 4 0) • Ar(fc). 
For the para-positronium (P) with Tp = we get 

sg } = ^Tr z 75 (l + 70)275(1 - To) = -2 
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kpky 1_ . 4k 2 



and 

Ep(fc) = ^ ■ ^Tr i 75 (l + 70)^75(1 - 7oK P = k 3*^ fc2 
For the ortho-positronium (V) with IV = jj we have 



1 

= 

and 



sf = ^Tr 7i (l + 70)7,(1 - 7o) = -2^ 



Sy(fc) = • -Tr 7,(1 + 7o)o- w ,7 i (l - 70)^ = 



fc2 4 - 1 ww/n- w»v ^3(k 2 + A; 2 ) 
The results are collected in Table 2. 

8 ^ 2 



3 k 2 + & 2 ' 

10 Appendix II 

Let us consider the contribution of the longitudinal part of the photon prop- 
agator to the integral (ITT]) . We have equality 

^(r)5 M /(z(r)) = |:/(s(r)) 

Then the term with the gauge dependent part d^d v D d looks like 

2 Si Sj 







e 





2 



<9ti<9t 2 



-D d (z®( 8i ) - z^\0)) + D d (z®(0) - z^(0)\ 

dk d(k 2 ) 1 - ^ x ~y) 



e 2 [D d (0) - D d {x - y)) = e 2 J 



(2^) 4 k 2 k 2 

This term does not contribute to the bound state mass and should be omitted. 
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11 Appendix III 

We consider the integral 

f [Dp -Jdr^p 2 (ryU(\p(r)\)]+^p] 

I = J da e = J -^-e 
p(o)=o, p(t)=o 

(41) 

where 

$[p] = -JJ dr x dr 2 A(\ Tl - r 2 |) • j&PW+Pte)) 
o 

The assumption is that the functional $ is small ||$[p]|| <C 1. We have in 
this case 

/ = J da = ATexp {-^ J da $ + 0($ 2 )| , N = J da. 
and the problem is to calculate the integral 

W = ($[p]> = ^jda<$> (42) 

1 /■ Dp -]dr[Ip 2 (r)-C/(|p(r)|)] 

p(o)=o, p(t)=o 
Let the Hamiltonian 

have the spectrum {E n } with the wave functions \& n (x) 

tf* n (x) = £ n * n (x). 
The time Green function can be represented in two forms 

f Dp -}dr[\p\r)~U{\p(r)\)] 



C 

P(t')=x',p(t)=x 

X;*n(x)e-(*-*'^*+(x'). (43) 
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The Green function satisfies the correlation for t > t" > t' 

GW(x,x') = J dy G t _ t »(x,y)G t «_ t /(y,x / ) 
We have for the function W 

t n 

W(t) = JdnJ dr 2 A(\n - t 2 \)H[t u t 2 ) 



o o 



with 



H(n,r 2 ) = ^jd Pl j dp 2 Gt-^pJjtPiGrL-^pJjSP'G^O) 

= E ^ II d Pl dp 2 e- E ^ x -^ n MKM)^ Pl 
. e -^(— )^ 2 ( Pl )^ 2 (p 2 )e^^e- E — ^ 3 ( p2 )^ 3 (0) 

_ J_ ^ e --B ni (X-T 1 )-£ n2 (ri-T 2 )-£; n 3T2 ^ ^ 
nin2".3 

• (| dp 1 *; i (/h)e i ' Pl *«(Pi)) (/ cZp 2 K> 2 )e^* n3 (p 2 )) *^(0) 

The function W(t) for £ — > oo behaves as W(t) ~ £Wo- It means that in the 
above stated sum in this limit only terms with n\ = n 3 = survive: 

H(n,T 2 ) ^Y, e ~ {En ~ Eo){Tl ~ T2) 



where 



d Pl n{ Pl )SP^ n { Pl )) (/ dp^M^^oM 
= E^ (E "- Eo)(T1 - T2) ^(^)^(-^ 

Con ^ = / rfx * (x)e^ x * n (x) 
Finally, we have for t — > oo 

W(t) = j dn jdr 2 A{\ n - r 2 |)Ee- (T1 -" )(B "- Eo) C 2 n (fj = t W 
oo n ^ ' 

Wo = Jdr A{r)^e-^-^C 0n C* n 
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Thus, 

/ = J da e" = ttjj(o) • e-KEo-Wo+o^)) 
11.1 Spherically symmetric potentials 

If the potential is spherically symmetric U = U(\p\) then the spectrum {E n ^} 
and the eigenfunctions are 

2£ + 1 

^nt m (p) = Rne(p)Y em (n), Yl Y im(n)Ytm(n) = — 

with ortho-normal conditions 

E *nUp)KtJfi) = Hp - p') = -Up ~ p')^ ~ n ') 

nlm P 

For radial functions we get (p = \p\) 

oo 

J dpp 2 R ne (p)R n , e (p) = 5 nn >, 
o 

OO 1 

J2R ne (p)R ne (pi) = -5(p-p>) 

n=0 P 

The form-factors looks like 

'k\ „ /k 



= J dpp 2 R 00 (p)R n i(p) J -JL • e^^ m (n) 



with 



— 



2 



y dpp 2 R o(p)Rni(p)je y 

^ ' 



tJ = \fe J '+* It 
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Finally we get 



W = Jdr A(r)£e^^*>(-1)<(2£+ 1)C& 



11.2 The Coulomb potential 

In the representation (|4ip the Hamiltonian is 



dp, 



1 

\p\ 



The Coulomb wave functions 

^nim{p) = Rne{p)Y em (n) 
are solutions of the non-relativistic Schrddinger equation (E n < 0) 
d 2 2d £{£+!) 1 



dp 2 ^ pdp p 2 
or according to [18] one can put p = 



H \E n i 

P 



Rne(p) = 
nr so that we have 



d 2 2d 
dr 2 r dr 



+ 1) n 1 



Kru(r) = 0, 



where the spectrum is 



n 



or E n t = E n 



An 2 



with n — 1, 2, 3, ... and £ = 1, ...,n — 1. 
Solutions are 



7t nl {p) 



1 



rr 



(2£+ 1)!\ 2(n-£ - 1)! 



e'^F l-n + £ + l,2£ + 2, 



C n i(k) = J dpp 2 Roo(p)R n i(p)je ( y 
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(4 + k 2 ) 2 
25QV2 k 2 



Several particular functions are 
Cio(fc) 16 

C 20 (k) 
C 30 (k) 
C 3 i(fc) 



(9 + 4P)3' ° 2lW 
432^3 k 2 (16 + 27k 2 ) 
(16 + 9/c 2 ) 4 

288^ k(16 + 27k 2 ) 
(16 + 9F) 4 ' 



128^ fc 
(9 + F) 3 



C 32 (£;) 



^6 6912 k 2 



5 (16 + 9k 2 Y 



11.3 "Time" potential 

We have the integral 



X 

e 4 

•Ma) = —I x (a), 
on 



Dp ~\ I ' drp 2 (T)+U[p,a} 



Ix(a) = J —e o 
U[p,a] = W[p;a]-W[p;0] 
In the paper [] it is shown that for small a 



(44) 



x 



U[p, a} = W [p; a\-W [p; 0] = ~ J dr\p(r)\. 

o 

The last integral corresponds to the one-dimensional non-relativistic quan- 
tum system with the Lagrangian 



L= P__ «i I 
L ' 4 8 IPI 



for which the Hamiltonian reads 



a , 



H = p z + -\p\ 



The Schrodinger equation looks as 

d 2 a , 



dp 2 8 



(p) = 5* (p) 



(45) 



(46) 
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Let us introduce 



2 

p = —, £ = ~r e 



then e is the eigenvalue of the equation 
d 2 



dv 2 



+ v 



Y(v,e) = eY(v,e), 



v e [o, oo). 



The non-normalized solution of the equation looks like 



Y(v,e) 



. V^eKi (l(v-e)i), 



v < e 



v > e 



The spectrum is defined by the equations 
d 



dv 



Y(v, e 2n ) 







for even states n — > 2n 



(47) 



v=0 



Y(0, e 2n +i) = for odd states n — > 2n + 1. 

The even eigenfunctions $ 2n (t>) = Y(v, e 2n ) should have n zeros and the odd 
eigenfunctions $ 2n+ i(t>) = Y(v, e 2n+ i) should have one zero for w = and n 
zeros for < v < oo. 

The wave functions are 



$n(p) 

The form-factors are defined like 

oo 

An (|) = J dp MP^n(p) 



p E (— oo, oo) 



i. 

«3 



1 

a a 



i 

Ct3 



dp J-OY — p,e y — p,e 



= < 



'OO 



/ 2 \ 

^ I dv cos (azqv) y ( v , e ) y (u, e 2n ) 



-7— L Jdt> sin lasqv) Y (v, e ) Y (v, e 2n+1 ) 
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with the symmetry condition 

A n (-asg) = {-l) n A n (a*q) , 
For «< 1 form-factors behave like 

A (aig) = 1 - at a , a = J ctoy Y = 0.374939, 

Ai (a-sq) = iJq a u a, = J dvv Y ^' = 0.862863, 

M (J g ) = aif a,, a 2 = / = 0.569709, 

~ V 2 

^ 3 = ^§g a 3 , a 3 = [dvv Y ^— ^/^J^ — = -0.0685378, 

The next integral for t — > oo behaves like 
W(t) 



t 71 

= JdnJ dr 2 H(|n-r 2 |) 
o o 

~ | d Pl | dp 2 G t _ T1 (0,p 1 )e i ^G ri _ T2 (p 1 ,p 2 )e^^G' T2 (p 2 ,0) 

t 71 

| dn / rfr 2 //(|n - r 2 |) ^ e -^-^)^- B °) A (a§g) A* K (-Jq) = t W 
oo K 





where 

oo 



>o = /dr ff(r)£e- r(VJ% ^ ^ (-ai?) 

o K 

Thus, for t — > oo one can write 



30 



References 

[1] Weinberg S. (1995) The Quantum Theory of Fields, Cambridge Univer- 
sity Press; 

[2] Wightman A.S. (1964) Introduction to some aspects of the relativis- 
tic dynamics of quantized fields, French Summer School of Theoretical 
Physics, Cargese, Corsica; 

[3] Wick G.C. (1954) Phys.Rev. 96:1124; Cutkosky R.E. (1954) Phys.Rev. 
96:1135; 

[4] N.Nakanishi, Suppl. of Pr.Theor.Phys., 43, 1 (1969); 

[5] G.V.Efimov, Few Body Syst, 33, 199 (2003); 

[6] Efimov G.V. (2007) Few Body Syst 41:157 (2007); 

[7] S.M.Dorkin, M.Beyer, S.S.Semikh and L.P.Kaptari, Few Body Sist, 42, 
1 (2008); 

[8] G.V.Efimov, hep-ph:08054337, (2008). 
[9] G.V.Efimov, Few Body Sist, 47, 137 (2010); 
[10] Dyson F.J. (1952) Phys. Rev. bf 85, 631; 

[11] W.E.Caswell and G.P.Lepage, Phys. Lett, bf 167B, 437 (1986); 

[12] De Alfaro V., Fubini S., Furlan G. and Rossetti C. (1973) Currents in 
Hadron Physics, American Elsevier Publishing Company, NY; 

[13] Ahiezer A.I. and Beresteckiy V.B. (1981) Quantum Electrodynamics, 
" Nauka" , Moscow; 

[14] Berestetskii V.B., Lifshitz E.M., Pitaevskii L.P. (1982) Quantum Elec- 
trodynamics, 2nd Edition, Pergamon Press, Oxford; 

[15] Sapirstein J.R., D.R.Yennie, (1999) Theory of Hydrogenic Bound States 
in Quantum Electrodynamics, Ed. T.Kinishita, Advanced Series on Di- 
rections in High Energy Physics, vol 7, World Scientific; 

[16] Eides M.I., Grotch H., Shelyuto V.A.: (2001) Physics Reports, 342:63; 



31 



[17] M.Dineykhan, G.V.Efimov, G.Ganbold and S.N.Nedelko, Oscillator 
Representation in Quantum Physics, m24, Springer- Verlag, Berlin Hei- 
delberg N.Y., 1995.; 

[18] L.D. Landau and E.L.Lifshic, Quantum mechanics, Moscow 



32 



Table I. QM and QFT. 



Quantum Mechanics 


Quantum Field Theory 


H = H + gH I 


(H + gH^E = Em E 
{^e} = {free particles} {bound states} 


H *<£> = E*f 
Fock space= {free particles} 
Hi is not operator on Fock space 


^(t) = e- iH ^- to ^(t ) 
Development in time 


S = \im t t ~^ +OD e iH ot e -iH(t-t ) e -iH t 
renormalization 
=>- S is operator on Fock space 
Relation between asymptotically free states 


intermediate particles ON mass shell 


intermediate particles OUT OF mass shell 


Relativistic corrections 


interaction is transmitted instantly 


retarded interaction 


Small corrections to H 
Elimination TIME out of Hamiltonian 
Effective theories 


Ain-+out(pi,P2,ki,k2) ~ M 2_( pi+p2 )2 

(0|J(x)J(0)|0> ~ e~ M W, \x\^oo 


Breit potential 
Nonrelativistic QED 


Bethe-Salpeter equation 


classification of states R 3 


classification of states R 4 = R 3 R] 
"Time" excitations =>■ abnormal states 
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Table 2. Quantum numbers of relativistic currents 



J 


{^Oji>) 


5 


L 


J 




J p 


s 


ft 1 \ 


(o-fe) 


1 


1 





+1 


0+ 




( T i \ 




-i 
l 


1 


(J 


+1 


IF 


V 


















(ifii'yij;) = 


=>- cr 


1 





1 


-1 


1 




(>7o7^) = 


=>- cr 


1 





1 


-1 


1 


1 


















(ll><Tijll>) => 


[cr X fe] 


1 


1 


1 


+1 


1+ 




(^757o^) 


1 











-1 





A 


















(^757^) =>• 


[cr X k] 


1 


1 


1 


+1 


1+ 


P 


(^75^) = 


1 











-1 
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Table 3. Functions Ep and E r (/c) 



J p 


r 


4°> 


Ar(fc) 


p(o-) 


«7 5 cos #p + 7570 sin 9 P 


-2 




v(i-) 


7j cos 6V + i7o7j sin 9 V 


-2% 


2 k 2 

3 k^+k'i 


5(0+) 


I cos 6* + 70 sin 9 










757 cos + % [7 x 7] sin 9 









Table 4. The function A n £ (a) 



(n£) 


(10) 


(20) 


(21) 


(30) 


(31) 


(32) 


A n , (0) 


1. 


0.0987 


0.0987 


0.0283 


0.0307 


0.00244 




0.9999 


0.09868 


0.09864 


0.02829 


0.03071 


0.002438 


A na (a) 


0.96707 


0.09387 


0.09617 


0.2683 


0.002986 


0.00238 


Anc(a) 


0.01421 


0.002089 


0.001173 


0.0006319 


0.0004049 


0.00002632 



Table 5. Norm N, 



n 




00 

N n = J dv & n (v) 






1.0188 


8.655 


1 


2.3381 


14.558 


2 


3.2482 


16.886 


3 


4.0879 


19.097 


4 


4.8201 


20.652 
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